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Abstract 

We discuss the recently discovered system SrCu2{BO^)2, a realization of an 
exactly solvable model proposed two decades earlier. We propose its interpretation 
as a Mott Hubbard insulator. The possible superconducting phase arising from 
doping is explored, and its nature as well as its importance for testing the RVB 
theory of superconductivity are discussed. 

Dedicated to Professor Bill Sutherland on occasion of his 6(f h birthday. 

1 Introduction 

Quantum spin systems are of great current interest, as shown by this symposium, with 
roots in two distinct sources. On the one hand the theory of model systems providing a 
rich variety of possibilities, and on the other, the field of synthetic materials, which has 
generated a vast number of systems, often close to theoretical models. As a result of this 
interplay, several interesting systems have been made in the laboratory, challenging our 
understanding by producing not only the expected, but also on occasion, the unexpected. 
Such a system that has caught attention recently is SrCu2(B0 3 )2, a two dimensional 
S = 1/2 isotropic Heisenberg antiferromagnet in two dimensions on a particular lattice 
with the property that it is solvable exactly for the ground state. Indeed it was solved 
two decades ago |1| by Sutherland and one of us. In this article we summarize the story 
so far, and also explore possible interesting physics that could arise if this system is 
doped. 

The situation of exactly solvable models in the area of statistical mechanics is rather 
limited. There is a general feeling that the special models are non-generic and rare, and 
hence somewhat ornamental. Enlightened opinion[§] has been more positive, and indeed 
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the role of some solvable models is very well recognized. In contrast, the situation in 
condensed matter physics is very positive. The interaction between new systems, new 
phenomena and novel concepts has been rewarding. The table below gives a few examples 
of popular systems, their realizations and the unique concepts associated with them. 



Year 


Model Systems 


Realization 


New Concepts 


1930 


one 

dimensional Heisen- 
berg Bethe AFM 


CPC, CuO chains 


Quantum Disorder, Spinons 


1968 


one dimensional Ma- 


(approximate 


Broken Discrete Symmetries 




jumdar Ghosh AFM 


mapping) CuGeO^ 


and Spinons 


1969- 


one dimensional 1/r 2 


Parametric Cor- 


Spinons, Unusual Statistics 


87 


Calogero, 

Sutherland, Haldane, 
Shastry systems 


relations in Quantum 
Chaos 




1969 


one dimensional 
Hubbard Model 


Benzene, Annulenes 


Spin Charge Separation, 
Holons, Spinons, SC Fluctu- 
ations from repulsion, Mott 
Hubbard Insulating state 


1988 


one dimen- 
sional Spin-1 Heisen- 
berg AFM, Affleck, 
Kennedy, Lieb and 
Tasaki chain 


Ni Chains, NENP 


Haldane Spin Gap 


1990 


one dimensional n leg 


Vanadates 


Integer vs non integer phe- 




Heisenberg Ladders 


CaV n 02n+l 


nomena, Superconductivity 
from doping Insulators 


1981 


two dimen- 
sional S=l/2 Shastry 
Sutherland model 


SrCu 2 (B0 3 ) 2 


Dimer states, Magnetization 
Plateaus... 



The 1/r 2 system in this table is different from the rest in that the physical realization 
comes from the world of quantum chaos, the continuum model is the description of 
parametric correlations in chaotic systems. The sole two dimensional system in the list 
is the main concern of this article. It has for long been unique in its very existence as 
a two dimensional member of the family of solvable models. It is particularly surprising 
since it is a model with essential simplicity as evidenced by the absence of crossed bonds. 
It is now even more remarkable in that nature finds a way of fulfilling the conditions 
for solvability in the compound SrCu2(B0 3 ) 2 - We discuss the origin of the model, its 
discovery in real life, some recent interesting developments in the physical properties, 
and some possible future directions. 
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1.1 Origin of the model 

In view of the enormous current interest in the problem, and also questions from col- 
leagues, it may not be inappropriate to say a few words on the Shastry Sutherland (SS) 
model on a special lattice, and how it came about. In 1980, I (BSS) joined the Univer- 
sity of Utah as a junior faculty member in the group consisting of Professors D C Mattis 
and B Sutherland. After an inspiring talk from Professor J R Schrieffer on polyacety- 
lene, I mentioned to Professor Sutherland that a clear magnetic analog of polyacetylene 
ground state is the Majumdar Ghosh (MG) model, the one dimensional Heisenberg with 
a second neighbour interaction half as strong as the first. Professor Majumdar, my PhD 
advisor at TIFR in Bombay earlier, had invented this model in an effort to go beyond 
the Bethe nearest neighbour antiferromagnet ( AFM). The model was well known to 
me, in spite of rather wise discouragement by Majumdar from working on Exactly Solv- 
able models, as a discipline unconnected with traditional topics in Solid State and Many 
Body Physics, in view of the almost zero probability of finding a new one! I remember 
being surprised that Sutherland, already then a sensei in the area of exactly solvable 
systems, had not come across this model! In that characteristically American way, there 
was no gap between learning of a new thing and getting excited and plunging into. I 
also caught the excitement that I had carried, but so far resisted (!). We first came up 
with the soliton excitations of the MG model, the so called spinons, as isolated unpaired 
spin 1/2 propagating objects in the midst of a sea of singlet dimerized spins. These 
were identical to the solitons of Schrieffer in spirit, but fractionalized the spin degrees of 
freedom rather than charge. Such excitations have since become a paradigm in the post 
High Tc language of strongly correlated systems, where the fixed singlets of MG give 
way to dancing singlets, the Resonating Valence Bond States envisioned by Anderson. 

In an effort to go beyond one dimension to higher dimensions, we tried various things. 
It was clear that a decomposition into triangles was the key to the MG model, and 
there was no essential reason why this had to be only one dimensional. The general 
point made was clear ||, the search for Integrable systems in higher dimensions is not 
very rewarding, the conditions for integrability seem hard to satisfy in higher than one 
dimension, however, the search for exactly solvable models ( for the ground state) is more 
promising a priori. In a d- dimensional Hilbert space there is a huge number ( ~ <i 2 /2)of, 
in general, non commuting operators that simultaneously share a given eigenstate, for 
example the dimer covering, and the search boils down to states and operators that 
satisfy the somewhat subjective criterion of "naturalness". As a result we pondered for 
several weeks on likely systems such as the two dimensional triangular lattice, where it 
became clear that no dimer like states work since the triangles share bonds with more 
than one other triangle. One needed a lattice where for a given triangle, no more than 
one bond is shared with another. This line of thinking led to the SS lattice shown in 
Fig.|(a). 

The proof of the ground state is simple and worth repeating if only briefly. The 
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Figure 1: The SS lattice. The angle 6 is the apical angle for triangles that are the 
building blocks of the lattice, and by continuously changing it, one generates different 
looking lattices with essentially identical topology, (a) represents the original choice of 
SS, and (c) the case closest to the Copper lattice of SrCu 2 (B0 3 )2- 



Hamiltonian can be written as a sum over triangles 

H = .m t E t = J I £ Si.Sj + 2a ss £ S*4J (1) 

\<i,j> <l,m> J 

where the subscript t refers to triangles, with H t = assSi-S2 + (5*1 + ^2). £3, 0:55 is the 
bond strength parameter, sites 1, 2 refer to the two sites on the diagonal and 3 the third 
site. Here and later we will denote the "dimer" bonds by I, m and the non dimer nearest 
neighbours as The first, and remarkable point is that the dimer state tp = Y\i <m [l, m ] 
is an eigenstate of H. Here the product runs over all dimers on the lattice, which must 
provide a covering of the lattice ( i.e. every lattice point must occur once and only once in 
the product). This happens because we can rearrange the operation of the Hamiltonian 
into two classes of terms, the wanted and the unanted terms. The wanted terms isolate 
the spin interactions on the dimer spins. Remarkably all unwanted terms have the form 
Sj.(Si + S < 2)[1,2], for appropriate indices, which vanishes on using the singlet property. 
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By Rayleigh Ritz variational principle Earner > Eq, but by the Anderson decomposition 
strategy we have a lower bound Eq > N t e t . Happily the upper and lower bounds coincide 
for ass — 1; an d we are guaranteed that this is the ground state. Later work|T5|] improves 
this lower bound on ass somewhat to ~ 0.74. 

These dimer ground states in two dimensions turn out to be surprisingly robust, 
for example the coupling constant ass is determined by an inequality rather than an 
equality, so there is an entire phase where the dimer states are the ground state, unlike 
the one dimensional case where one has a solution at isolated points only. This clearly 
greatly increases the probability of finding such models realized in nature, whereas in one 
dimension we should only expect proximities. Further the ground state is insensitive to 
the spin space isotropy of the underlying Hamiltonian, and one has the strange situation 
where the ground state has a greater symmetry than of the Hamiltonian! The phrase 
"superstability" || describes this kind of robustness shared by most of the dimer ground 
state systems. An example of robustness comes from later in the story, where we find 
that the ground state of stacks of the SS lattice, rotated by it/ 2 and coupled by vertical 
spin interactions, a model that describes the real 3-d material SrCu2(BOs)2, "magically" 
manages to have the dimer state as the true ground state[|J! 

1.2 The system SrCu 2 (B0 3 ) 2 

Almost two decades later, Kageyama and coworkers at the ISSP in Tokyo found that 
SrCu2{BO^)2 1 synthesized earlier in 1991 by Smith and Keszler, had very unusual prop- 
erties. The spin 1/2 moments of Copper living in well isolated two dimensional layers 
seemed to lock up into singlets, and a clear spin gap behaviour was observed by NMR. 
They concluded that this is a unique system, the first truly two dimensional spin gapped 
system with S = 1/2. The data was analyzed by Mihayara and Ueda (MU)[|J, who 
realized that the physical system was describable by an exactly solvable model. They 
proposed the model, found its solution, and then realized that it was essentially (topo- 
logically) the same model as SS, but looked different due to the details of the lattice. 
The Copper lattice is shown in Fig.[TJ(c), and by opening up the angle 9 continuously, one 
reaches the SS lattice ( with 9 = ir/2 ). An intermediate value of 9 = 7r/4 in Fig.|l](b) 
aids the imagination. Changing the angle 9 clearly preserves the orthogonality of the 
"dimer" bonds but changes the bond length relative to the inter dimer bond lengths. This 
is crucial, since the criterion for solvability ass >~ 0.74 becomes realizable only in this 
picture. A nice visualization of this deformation is available courtesy Dr H Kageyama 
at |http:/ /www. issp.u-tokyo. ac.jp/labs/mdcl/ueda/kage/head. html| . The Hamiltonian is 



written by MU and some recent papers as H = J[aMuJ2i<j Si-Sj + J2i< m Si-S m ], with 
J' = Jauu-, an d hence we clearly have ass • °-mu — §■ The current estimates of 
(J, auu) using experimental data on SrCu2{BO^)2 range from ||] (85°K, 0.635) to || 
(71.5°/^, 0.603). Thus ass ~ 0.78 is perilously close to, but on the safe side of the phase 
boundary at ~ 0.74. 
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1.3 Recent Developments 

We next mention a few of the very large number of papers that have been written 
recently, with apologies in advance for possible incompleteness. After the discovery 
of the spin gap, neutron scattering has confirmed the absence of magnetic LRO and 
inelastic scattering has given clear indication of a flat dispersionless triplet excitation 
mode at about 3 meV, as well as of many branches of dispersing bound states of triplet 
excitations || . NMR experiments were the first to show the spin gap|7j at about 30 K. 
ESR experiments show the presence of a second gap at about 4.7 meV which implies 
a substantial binding energy of two triplets ||. Raman studies show a singlet bound 
state at about 3.7 meV |J. The magnetic exchange constants, as mentioned are in the 
60-80 K range. This is convenient for exploring with available pulsed high magnetic 
field experiments, which reveal |7|, |12| the surprising existence of magnetic plateaus at 



M/M s = 1/4, 1/3, 1/8... There is interesting data on the effect of magnetic excitations 



on phononic thermal conductivity! 10 



Thus a large set of experiments have already been done, and provide many constraints 
on the theory. We should mention that the knowledge of the exact solution of the ground 
state does not give much insight into the excitations in this class of systems. One knows 
that in general terms, the singlet dimers can be broken into triplets, and that isolated 
triplets find it hard to propagate on the SS lattice due to its topology. This leads to 
flat bands of triplet excitations, i.e. very massive objects, consistent with neutron data. 
Pairs of triplets, however, escape the topological constraints, much as holes in the Neel 
Antiferromagnet, and move about quite freely. Thus one has kinetic binding and the 
bound state has substantial dispersion [p4j|. 

In a sense the unexpected and new physics so far has been the presence of these 
plateaus [Q. These are unique in that they are the first two dimensional plateaus seen, 
and have attracted considerable interest. Several possible scenarios have been suggested 
to explain these. One picture is that of massive triplet excitations acting as hardcore 
bosons, that hop as well as interact. The effective interactions are strong due to large 
mass, and Wigner crystallization is proposed to explain the plateaus |13|. Alternatively 
one can view this plateau formation as the Quantum Hall Effect of hard core bosons, 
and a Chern Simons type field theory provides a fair description. The structure of the 
Hofstadter spectrum on the SS lattice is reflected in the plateaus [14]]. At the moment 
it is not easy to reach a conclusion as to the best interpretation, especially since the 
experiments do not show plateaus that have anything like the precision of the Quantum 
Hall Effect. 

The phase diagram at zero field has come in for close scrutiny by several authors, 
using series expansion ideas (15], [16|], rigorous bounds |17|], large N field theories fl8|], as 



well as effective field theories of bosonic dimers ]|T9fl . Exotic and unusual intermediate 
phases are suggested by these studies, including a "topologically ordered phase with 
deconfined 5 = 1/2 spinons, which should give rise to an exotic superconductor with 
anomalous flux properties under doping". An early paper by Albrecht and Mila|20| 
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Figure 2: The phase diagram of Ref. in the large N limit. The abscissa may be read 
as 2a ss /(l + 2a ss )- 



discusses the transition between the dimer and AFM phases using Schwinger Bosons, 
and concludes that it should be first order. A qualitatively new intermediate spin liquid 
phase with power law correlations has been proposed by Koga and Kawakami||T5|. The 
phase diagram in Fig.|2| is from reference for a large N theory, and represents a 
possible set of phases with various kinds of magnetic order. The Koga-Kawakami phase 
may be viewed as the SRO (ir, q) phase. 

We should also mention new theoretical models that are generalizations of the SS 



ideas to higher dimensions and other systems [221, BTJ. The beautiful model of M tiller 



Hartmann, Singh, Knetter and Uhrig [22| has a new set of exchange interactions added to 



the SS model, and had a very rich set of constants of motion. It is very tractable, giving 
rise to magnetization plateaus that are similar to but not identical to the experimental 
ones. 
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2 Mean Field Theory of the Doped Dimer Super- 
conductor, a Possible test of RVB 

In this section we discuss the possibility of doping the dimer state, and what one 
may expect from it. Firstly we remark that the insulating dimerized ground state of 
SrCu2(BOs)2 may usefully be considered as a Mott Insulating state of an underlying 
Hubbard model. To see this consider the Hubbard model on the SS lattice with 

H = -t J2 ( c L c jv + h - c -) -a f E ( c l,* c m,« + h - c -) + U H n r,~\n rA (2) 

<ij>,a <l,m>,a r 

where the (somewhat overused) symbol a represents the ratio of hoppings on the two 
kinds of bonds, and the Hubbard interaction term sums over all types of sites. Clearly 
the superexchange argument fixes it in terms of the ratio of the exchange parameters 
via a 2 = 2ass, and we note that a = ±1.25 using the insulating state estimates, with 
the sign undetermined. In the non-interacting limit the band structure is interesting, 
we have four subbands, with the extrema of two of them touching quadratically at the 
zone center. At half filling, one has four electrons per unit cell and the system is a 
semi-metal with a finite density of states, and thus it has typical metallic behaviour such 
as a linear specific heat. A parallel may be drawn with the semi-metallicity of graphite 
on the hexagonal lattice with two electrons per unit cell and also of a fiduciary MgB 2 
with well separated planes. In the case of graphite however, one has a "Dirac like" linear 
spectrum, and hence the density of states near the "fermi point" , i.e the contact point 
vanishes. 

This semi-metal becomes an insulator at large enough U, undergoing a transition to 
the dimerized state that does not break rotation invariance nor the lattice translation 
symmetry, and may be called a Mott transition in the same sense as that of the one 
dimensional Hubbard model at half filling at infinitesimal U. Since the large U behaviour 
is exactly known, namely the dimer ground state, further terms in the t/U expansion 
beyond superexchange should be useful in elucidating the nature of the Mott transition 



here |23|| . This transition has not yet been studied in literature. Starting from the 
semi-metal and turning on U, one may either have a level crossing transition to the 
gapped insulator, or more interestingly a continuous opening of the charge gap. In the 
gapped insulating phase, the four spin correlation function pertaining to dimer order 
< SaSbScSd >~< S a Sb >< S c Sd >, and thus there is ODLRO in this correlation 
function without any obvious symmetry that is broken. With this distinction, without 
necessarily a major difference, we may refer to SrCu2{BO%)2 as a Mott insulator. 

Having this realization of the Mott Hubbard insulator, we naturally enquire if the 
philosophy of the RVB theory of Superconductivity due to Anderson applies here. This 
theory is built upon the idea that repulsive interactions of the Hubbard type lead to 
superconductivity via the intermediate step of superexchange, or Heisenberg interactions 
in the insulating state. The superexchange leads to singlet pairing between electrons of 
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opposite spin, and these pairs are analogous to the Cooper pairs, but are localized due to 
the Mott-Hubbard gap. Under doping the Mott-Hubbard gap collapses, these preexisting 
pairs propagate freely, and lead to superconductivity. In the present case, the Cooper 
pairs at half filling should be viewed as the dimer- singlets, which on doping should move 
around by the same logic, and lead to superconductivity. Since the values of exchange 
are smaller by an order of magnitude from those in High Tc systems, we expect lower 
Tc's, say tens of degrees K, but accompanied by the characteristic signature of singlet 
pairing and also of definite phase relations of Cooper pairs on the bonds, analogous to the 
d-wave pairing. While this theory is remarkably effective in providing a comprehensive 
view point, it still lacks unambiguous experimental support or a rigorous mathematical 
foundation, and one would welcome other supports to its validity or otherwise. In this 
context we work out in this section the mean field theory of a fiducial doped SrCu2(BO)s, 
and calculate some characteristics of the proposed superconducting compound. 

Before doing so, let us note that doping can be of either chemical type, as in say 
Sri_ x M x Cu2(B03)2 with a monovalent alkali Mora trivalent lanthanide. However, one 
interesting possibility is suggested by the comparison of MgB2 with graphite. One learns 
that MgB2 is isoelectronic with graphite, but avoids being a semimetal by dispersing 
the bands in the direction transverse to the two dimensional sheets, it self-dopes by 
decreasing the transverse lattice constant. It is possible that a divalent element like Mg 
in place of Sr with a smaller ionic radius could decrease the transverse lattice constant 
of SrCu2{BOzj2 sufficiently so that it would have substantial transverse dispersion. 
We should clarify that unlike MgBi which appears to be a case of phonon mediated 
superconductivity ||24|| , we are examining the case for a non phononic mechanism, the 



RVB mechanism for doped SrGu^iBO^- Gate charging might be another attractive 
possibility. We now turn to the calculation proper. 

2.1 RVB type mean field theory on the SS lattice 

We next present the mean field theory of a t — J type model on the SS lattice. The 
hopping amplitudes and the exchange integrals on SS lattice are as shown in Fig.(^). 
The nearest neighbour (n.n.) hopping amplitude is — t, and the next nearest neighbour 
(n.n.n.) hopping amplitude is —at where a is a dimensionless number. The exchange 
couplings, J and J' along the n.n. and n.n.n. directions respectively, are such that 
J' = a 2 J, as governed by the large U physics of the Hubbard model on SS lattice. 
Thus a 2 = 2ass of the previous section. The t — J type model, thus arrived at is an 
appropriate generalization of the original SS model, in order to deal with doping. In 
the following, we will first describe the tight binding band structure of the free electrons 
on the SS lattice. Then, we will do the mean field theory of the interacting model and 
discuss its implications for superconductivity, in a manner analogous to the early RVB 
mean field theories of t — J model on a square lattice done in the context of high-Tp 
superconductivity^, EH . 
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Figure 3: The Shastry-Sutherland lattice. Also shown is the labeling of sites in a unit 
cell of SS lattice, as used in the text. 



The band-structure : The SS lattice has a periodicity of 2a a , both along x as well as y 
directions, where a Q is the lattice constant. With each site contributing just one relevant 
orbital, the tight-binding model on SS lattice is described by a four band hamiltonian 
given below. 

c (k) 



H t = J2[ct(k) d(k) c|(k) 4(k)l T(k) 



k,<T 



C1 (k) 

c 2 (k) 
c 3 (k) 



(3) 



Here, a =| or |, and the wave-vector, k = (k x , k y ), is such that ^ < k x , k y < The 
subscripts, 0, 1, 2, 3, refer to four different site within a unit cell. The dispersion matrix, 
T(k), is a 4 x 4 hermitian matrix as given below. 



T(k) = -t 



2 cos(k x a ) 2 cos(k y a ) a e l( - kx ky ^ a ° 

2 cos(k x a Q ) a e i( - kx+k ^ a ° 2 cos(k y a ) 



2 cos(k y a ) a e~^ kx+k ^ a ° 



a e 



-i(k x k y )a 2 cos(k y a ) 2 cos(k x a ) 



2 cos(k x a ) 







(4) 



The band-structure for \a\ = 1.25 is shown in Fig.(|J). This value of a is taken from 
the studies on SrCu2(B03)2, where the values of J and J' are extracted by fitting the 
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experimental data with the orthogonal dimer model. What is known from the experi- 
ments is a 2 , and not a. This leaves us with the ambiguity of sign of a, and hence we 
have considered both positive as well negative values of a. 



a = -1.25 




.4 t 1 1 ] 

1(0,0) M (It/2, 0) X (t/2, Jt/2) no,0) 



a = 1.25 




1(0,0) M (Jt/2,0) X (71/2, 71/2) T(0,0) 



Figure 4: The band-structure of free electrons on the SS lattice. The wave vectors are 
written in units of — . Notice that the band-structure is odd with respect to a. 

Let us make a few essential observations regarding the band-structure. Firstly, the 
system is a semi-metal at half filling, since the middle two bands touch each other at the 
zone centre. Secondly, there is a band which is flat along the Xr symmetry direction in 
the Brillouin Zone. This band gives rise to a severe van Hove singularity at a. Thirdly, 
the values of band energies at zone centre are (—4 — a), a, a and (4 — a). For \a\ > 2, the 
middle two bands no more touch each other, and there is a finite band gap which makes 
it a band insulator at half filling. Since a for the material of real interest is roughly 1.25, 
we have not tried to discuss other values of a. 

Fig.(|5|) shows the non-interacting single particle density of states on SS lattice for 
both negative as well positive values of a. When we hole-dope the system to take it 
away from half filling, it is expected to behave differently for positive and negative a, 
since the flat band influences the case only when a is negative. 

The mean field hamiltonian : The t — J hamiltonian on SS lattice can be written 
as: 

H = VH t V + ffj - t/i E E 4,M c pAV (5) 

k,o- p=0 

The first term in H accounts for the projected hopping. It is essentially H t as given in 
equation^), but with projection operator V which suppresses the double occupancy of 
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a = -1.25 a =1.25 




E (in units of t) E (in units of t) 



Figure 5: The density of single-particle electronic states on Shastry- Sutherland lattice. 

any site (due to large Hubbard U). At a simple level, the effect of V can be brought in by 
replacing t by St. Here, S(= 1 — n), is the number of holes per site, and n is the electron 
filling per site. The last term is the chemical potential, fit, times the total number of 
electrons. Here, p is the site (or the orbital) label within a unit cell. The second term in 
equation^), Hj, which accounts for the interaction among electrons can be written as: 

^ = 4E+« 2 E}fs(r).S(rO-^>) (6) 

(n.n. n.n.n.) \ / 

Here, r and r' are the site labels, and h(r) denotes the number operator at site r. 
The summation is pairwise in r, r'. The primed summation denotes the sum of only 
those pairs of n.n.n. sites which are allowed by the connectivity of the SS lattice. The 
operator, (S(r) ■ S(r') — n(r)n(r')/4), can also be written as — |bt(r, r')b(r, r'), which 
provides the basis for mean field decoupling of Hj in the off-diagonal channel. The 
operator, b(r,r') = C|_(r)cj(r') — cj (r)cj(r'), is the singlet bond operator. 

Let us define an off-diagonal or the pairing mean field, (b(r, r')), in the following way. 

{Ae i9x for r — r' = ±ax 

Ae ^ for r - r' = ±ay 

A ' e ie x+y for r - r' = ±a(x + y) 1 1 

A'e^-y for r — r' = ±a(x — y) 

The phases, 9 X , 9 y , 9 x+y and 9 x - y , as well as the amplitudes, A and A', are all independent 
of the coordinates. Hence, we are considering a uniform case. With this choice of the 
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order parameter, we decouple Hj. The corresponding mean field hamiltonian can be 
written as: 

H MF = H t + Hf F + L(4JA 2 + J' A' 2 ) (8) 

where L is the number of unit cells. In order to write H t and Hf F conveniently, we 
introduce a notation. Let us define the Nambu operators, \I/ C f (k) and \E^(— k) in the 
following way. 



# CT (k) 



-k) 



cot(k) 

c it( k ) 
c 2T (k) 

C3t(k) 



[coi(-k) c u (-k) c 2i (-k) c 3i (-k)] 



(9) 



(10) 



The subscripts, C and R, indicate that ^c^k) is a column vector and ^^(— k) is a row 
vector. In this notation, H t can be written as : 



H t = Y, tr{T(-k)}+[^ T (k) tf^-k) 



T(k) 
-T(k) 



*ct(k) 



(11) 



Here, T(k) is essentially same as the dispersion matrix, T(k), except that the chemical 
potential forms its diagonal elements, and all the off-diagonal entries have a factor of 
hole doping, 5, in order to account for the projection. 



T(k) 



fx 25 cos(k x a ) 25 cos(k y a ) 5a e i( ~ kx ~ ky ^ ao 

25 cos(k x a ) fx 5a e i( > kx+k ^ a ° 25 cos(k y a ) 

25 cos(k y a ) 5a e ~ i( - kx+ky ^ a ° fx 25 cos(k x a ) 

5a e~ i( ~ kx ~ k ^ a ° 25 cos(k y a ) 25 cos(k x a ) fx 



With the same notation, Hf F can be written as 



jjMF _ 



*kr(k) ^(-k) 



D(k) 
Dtfk) 



^ CT (k) 
-k) 



Hi* 



where D(k) is a non-hermitian 4x4 matrix as given below. 



D(k) = 



JAe ie *cos(k x a ) 

JAe ie *cos(k x a ) J^i e ie x+VQ i{k x +k v )a jAe i9y cos{k y a ) 

JAe if) ycos(k y a ) £AL e ie x+ye -i(k x +k y )a q JAe i9x cos(k x a ) 
£AL e i0 x -y e -i(k x -k y )a o jAe Wy cos(k y a ) JAe i9x cos(k x a ) 

(14) 



JAe i9y cos(k y a ) 



J' A' 



(12) 



(13) 



■J.0 x —y pi(&x ky) 
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Here, J' = a 2 J as mentioned earlier. Finally, we write the H MF as : 



H MF = E[^ T (k) <M-k) 



T(k) D(k) 
Dt(k) -f(k) 



-k) 



+L (4 JA 2 + J'A' 2 - 4t//) 



(15) 



Let us denote the matrix 



by A(k). It is an 8 x 8 symplectic, hermitian 



T(k) D(k) 
Dt(k) -T(k) 

matrix whose eigenvalues are real and occur in pairs. That is, an eigenvalue's negative 
is also an eigenvalue. 



The mean field free energy and the self-consistent equations : The grand 
canonical free energy, $, at a given temperature T, for the mean field hamiltonian 
described above is, 



$ = AL JA 2 + 



J'A 



i2 



^) " EE {^(k) + |Log (1 + e"^») } (16) 



Here, (3 = l/k B T, and |£^(k),j = 1,4} are the positive eigenvalues of A(k). Let us 
put t = 1. Now, all the energies (or parameters with units of energy) are in the units of 
t. We find the self-consistent equations for A and A' by minimizing $ with respect to 
A and A'. These are as follows. 



A 
A' 



--EE^tanh^' 
J'iL^U dA' \ 2 j 



(17) 



Since d$>/dfi = —N, where N is the total number of electrons, we get the following 
equation for the chemical potential. 



8 



A T E E 

k j 



dE+(k) 



4L ^ dfj, 



tanh 



(19) 



The hole doping, 5 = 1 — N/4L. Solving these sets of equations self-consistently gives us 
A, A' and \x as a function of 5, for given values of a, J, (3 and the phase angles 9 X etc. 



The results of the mean field theory : We solve equations (IT?]), (|18D and fll9D 
self-consistently for different values of 5. We are interested in both the hole as well as 
electron doping for a given a. It is clear from the band structure that the hole doping 
for a is same as the electron doping for —a. Therefore, we have considered only the hole 
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doping for both positive as well as negative a. In all our computations, we use t = 1, 
J = 0.3 and |a| = 1.25. The value of J for SrCu2(BOs)2 is roughly 70°K. The ratio of J 
to t is tentative, and taken to be roughly same as that for the high-T(7 superconductors. 
Though we have four phases, 9 X , 9 y , 9 x+y and 6 x _ y , only three relative phases are relevant. 
Therefore, we keep 9 X = 0. We have to find out those values of 9 y etc. for which the free 
energy is minimized, and see how things evolve as a function of 5. 

Let us discuss the zero temperature {(3 — > oo) case first. Fig. (P) shows the variation of 
A and A' with respect to 5, for a = —1.25 and 1.25 at zero temperature. For a = —1.25, 
the minimum of free energy occurs for 9 x+y — 6 x - y = vr regardless of the values of 9 X and 
9 y , and A is identically zero. For a = 1.25, the minimum of free energy corresponds to 
9 X = 0, 9 y = 7r, 9 x+y = and 9 x - y = tt. It is a weak minimum as many other choices 
of the phases have similar values of the free energy. Nevertheless, this choice of phases 
appears to be the minimum. It is interesting to note that for 5 = 0, the diagonal bond 
order parameter, A', is one and A is zero (and is independent of the phases 9 X , 9 y , 9 x+y 
and 9 x _ y \ Thus, the RVB mean field theory at half filling exactly reproduces the known 
dimer ground state of the SS model. 



a = -1.25 a =1.25 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



5 5 



Figure 6: The variation of A and A' with doping, S, in the ground state. For a = —1.25, 
A is not shown, since it is identically zero. 

To consider superconductivity in our mean field theory, we define a physical order 
parameter, Asc = -FbAmf- Here, Amf is the mean-field order parameter (A or A' 
whichever is larger for a given doping), and Fb is a bosonic mean field. Such an order 
parameter can be understood in the framework of slave boson approach p7|, where the off- 
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diagonal order parameter of the physical electrons is described as \bibjf} a fj a ,j. Here 6's 
are the slave boson fields, and /'s are the fermionic objects. In a mean field decoupled 
theory, this is like (bibj) {^fi a fj a i^ = -FbA M f- The bosonic order parameter, F B} is a 
function of temperature and doping, and goes roughly like 5. The superconducting 
transition temperature, Tsc, is the temperature where either Amf or Fb vanishes first 
while increasing the temperature. For low doping, Amf is large, therefore, the Tsc is 
same as the bose condensation temperature, T BC , for the bosonic field. Some estimates 
of Tbc have been made earlier while studying t — J model in the context of the high-T^ 



superconductivity p7 1 . We roughly estimate it by considering an approximate dispersion 



of the form, k\ + k\ + -k 2 z , with the z-axis anisotropy 7 ~ 30. We get T B c ~ 47rp*5(l — 
5)/[2+Log(47/7r)]. Here, p* is the density of states at the energy where two middle bands 
touch, from the side where dispersion is quadratic, and is a measure of the curvature of 
the band. For \a\ = 1.25, p* « 0.1, thus T BC ~ 0.225(1 - 5). 

One comment should be made regarding the interpretation of the result A = and 
A' 7^ (for a < 0). While at half filling this implied the dimerized insulating state, away 
from half filling it must be interpreted as superconductivity. The BCS type wavefunction 
implies the fermion pairing in real space, 



t. „t \ _ v-^k.^-r,). 



y „ 

k y/Al + (e k - /i) 2 

and it extends over a range of lattice constants (due to the non-trivial k dependence of 
tk away from 5 = 0), despite the mean field Hamiltonian having n.n. pairing only. A 
similar remark holds for the four fermi operator that determines the superconducting 

ODLRO of Yang, namely (cjt c ii c -?'-l- Ci '1') ^ ® ^ or l r * ~~ r il >> ^ - 

Fig.(|^) shows the phase diagram in the T-5 plane, as estimated from our RVB mean 
field theory. The temperature, T M f, is where A M f vanishes. The estimated T B c, and 
the computed Tmf are plotted as a function of 5. The common region under these two 
curves is the superconducting phase bounded by critical lines. As usual all the remining 
lines should be viewed as crossover lines rather than critical lines. Among the remaining 
three regions of T-5 phase diagram, the low doping region below T MF and above T BC 
is the spin gap phase with a suppressed density of states manifested in susceptibility as 
well as optical conductivity. Similarly, the high doping region is the normal fermi liquid. 
There is a region which is usually referred to as the strange metal phase, as shown in 
the Fig.(0) with linear resistivity. Also, the phase diagram is similar for both positive 
as well negative values of a. From this phase diagram, we estimate the optimal value of 
the superconducting transition temperature, T c ~ 10°K. 

In conclusion, the system considered here has a rather rich history. It may also have 
an important future since under doping it might be the much sought after low T c RVB 
superconductor, with linear resistivity down to 10 K and other such exotic properties, 
rather than a conventional phononic BCS superconductor. 
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a = -1.25 



a = 1.25 



0.2 



0.15 



0.05 



0.2 




Figure 7: The phase diagram for the negative as well positive values of a. The lines 
of the estimated Bose condensation temperature, Tbc, and the computed mean field 
temperature, Tmf, divide the T-S plane into four physically distinct regions. Each of 
these is appropriately labeled. 
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